CMB polarization map derived from the WMAP 5 year data through Harmonic 

Internal Linear Combination 
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We have derived whole-sky CMB polarization maps from the WMAP 5 year polarization data, 
using the Harmonic Internal Linear Combination (HILC) method. Our HILC method incorporates 
spatial variability of linear weights in a natural way and yields continuous linear weights over the 
entire sky. To estimate the power spectrum of HILC maps, we have derived a unbiased quadratic 
estimator, which is similar to the WMAP team's cross power estimator, but in a more convenient 
form for HILC maps. From our CMB polarization map, we have obtained TE correlation and E 
mode power spectra without applying any mask. They are similar to the WMAP team's estimation 
and consistent with the WMAP best-fit ACDM model. Foreground reduction by HILC method is 
more effective for high resolution and low noise data. Hence, our HILC method will enable effective 
foreground reduction in polarization data from the Planck surveyor. 

PACS numbers; 98.70.Vc, 98.80.Es 



I. INTRODUCTION 

The Cosmic Microwave Background (CMB) is ex- 
pected to be linearly polarized by Thompson scatter- 
ing at the last scattering surface and after re-ionization. 
Most useful information on re-ionization and primordial 
gravitational wave can be obtained from CMB polariza- 
tion on large angular scales P, 0, H, 01 ■ 

Foregrounds degrade the attainable accuracy of cosmo- 
logical information Q and Galactic foregrounds are par- 
ticularly significant on large angular scales [1, 0]- Hence, 
the ability to clean foreground contamination without 
relying on masking is of the utmost importance for the 
study of CMB polarization on large angular scales. 

Recently, multi-frequency polarization data for the 
whole sky has become available from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) [1, Q. For its cos- 
mological analysis, the WMAP team used foreground- 
reduced maps obtained by a template fitting method. 
Due to heavy foreground contamination within their 
Galactic mask, the WMAP team's foreground-reduced 
polarization maps are not suitable as whole-sky CMB 
polarization maps. The WMAP team has also pro- 
duced a low resolution CMB polarization map by Markov 
Chain Monte Carlo (MCMC) method ^9^. This MCMC 
map contains heavy foreground contamination within the 
Galactic cut @, making it unsuitable for whole-sky po- 
larization maps. In contrast, the Internal Linear Com- 
bination (ILC) method is one of best blind approaches 
available in spite of a few drawbacks. As a result, several 
variants of ILC methods have been developed and im- 
plemented to construct a whole-sky CMB temperature 
map [ll[ll|,[ll[ll[ll[ll. Meanwhile, whole-sky CMB 
polarization maps are quite scarce at the moment. We 
have derived a whole-sky CMB polarization map from the 
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WMAP 5 year polarization data through the Harmonic 
Internal Linear Combination (HILC) method, which is a 
natural extension of our previous efforts to derive a CMB 
temperature map [l3[. 

The outline of this paper is as follows. In Sec. [TTl we 
discuss the Stokes parameters in the context of an all-sky 
analysis. In Sec. Illli we briefiy discuss foreground reduc- 
tion method with multi-frequency maps. In Sec. IIVI we 
derive equations and solutions. In Sec. IVIl we present 
the results of applying on the WMAP five year polariza- 
tion data. In Sec. IVIIl we briefly discuss a few other 
foreground reduction methods and compare them. In 
Section rVIIIi we summarize our investigation and present 
our conclusions. In Appendix [Xl we discuss the unbiased 
quadratic estimator of power spectrum that we shall be 
using, which is similar to the WMAP team's cross power 
estimator, and in more convenient form for an HILC im- 
plementation. 



II. STOKES PARAMETERS 

The Stokes parameters describe the state of polariza- 
tion 

[1, 

\vf\ . For an all-sky analysis, they are measured 
in reference to (ee,e^) [H, 0]: 

Q = {El -El), 
U = {2EgE^), 

where (. . .) indicates average over a period of the elec- 
tromagnetic waves, and eg and are basis vectors of 
spherical coordinates. 

All-sky Stokes parameters are decomposed into spin 
±2 spherical harmonics 2]: 

Q{n)±iU{n) ^ ^ a±2,im ±2Yim(n). (1) 

l.ni 

The spin ±2 spherical harmonic coefficients a±2.im a-re 
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further decomposed into E and B mode [3, [T^ : 

a±2,hn = —{aE,im^ia,B,im), (2) 

where 

aE.lm — \-^) 0-E,l-nn 
aB,lm = (-l)™aB,i-m- 



III. FOREGROUND REDUCTION WITH 
MULTI-FREQUENCY MAPS 

Polarization at a frequency Vk and pixel x is as follows: 
(3(x, ^-fc) ±iU{-x,iyk) = 

(3cmb(x) ± «J7cmb(x) + Qfg(x, J/fc) ± iUlg{-K, Vk) 

(x, Vk) ± «J7„oiso(x, I^fc), (3) 

where 'fg' denotes the composite foreground signal. A 
natural candidate for the estimator of the CMB polar- 
ization map is a linear combination of multi-frequency 
maps, which is as follows: 

^u.^-(x) (g(x,i.fc)±i[/(x,z.fc)). 



and hence the linear combination map of minimum vari- 
ance contains minimum residual foregrounds. However, 
the variance minimization on real-world data proceeds 
in the way to maximize the cancellation between the 
residual foreground and CMB [H, [iB|, since the cross 
terms (Re [((3cinb(x) ± i[/cmb(x)) C(x)]) completely van- 
ishes only in the limit of infinite number of pixels. The 
method to reduce the cross term effect is discussed in 

The frequency spectrum of real-world foregrounds 
varies with their positions (see [l^ for a recent treat- 
ment), and a spectral index variation as small as ~ 10% 
can have a substantial impact on the linear combina- 
tion of multi-frequency maps 20]. Therefore, the linear 
weights should possess some spatial variability. At the 
same time, the linear weights of minimum foregrounds 
are expected to be spatially coherent on small angular 
scales to a good approximation. Hence, we assume the 
w^{6,(j)) to be some angular functions, which contain 
spherical harmonics of multipoles up to some cutoff mul- 
tipole ^cutoff- The linear combination map formed from 
multi-frequency maps 

Q{0,<j))±iU{e,^) 

= J2 ^'(^^ ^) ^k) ± tuie, </., i^k)) , (7) 



To keep the CMB unchanged, we choose linear weights 
such that the sum of linear weights over frequency chan- 
nels is equal to unity: 



(4) 



With Eq. [3] and 21 it is straightforward to show that 
E«;'=(x) (Q(x,;/fe)±z[/(x,i.fc)) = 

k 

(3cmb(x) ± iC/cmb(x) + ((x), (5) 

where we have neglected noise and defined C(x) by 

C(x) = XI "'''W (gfg(x,j/fc) ±ic/fg(x,i/fc)) . 
fe 

The variance of the linear combination map is 

Y.'^^i^) (Q(x,l^fc)±l[/(x,l.fc)) 

fc 

« C2 + 2 (Re [(Qcmb(x) ± z?7c„,b(x)) C(x)]) 
+ (IC(x)P), (6) 

where (. . .) denotes average over the whole sky. The 
term in Eq. [6] is the variance of CMB polarization 
signal and is, therefore, independent of linear weights 
because of Eq. [51 Since there is no correlation among 
CMB and foregrounds, we find a"^ ^ + (|C(x)P), 



can be rewritten in spherical harmonics through the 
Clebsch-Gordon relation: 



a±2,LM 

(-1) 



(8) 
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±2 =f2 J / ^ "^'™ "±2,i'm'; 



where 



a±2.LM = / ±2YIm{o, 4>) {Q{o, 0) ± iu{e, </))) dil, 



2V„Ao, 0) {Q{e, cj), Vk) ± iu{e, 0, vk)) dn. 



'■±2,l'm' — / ±2 



The variance of the linear combination map is equiva- 
lently given by 



LM 

where L is bounded by the triangular inequalities: 

\l-l'\ <= L <=l + l'. (9) 

The constraint J^i'^'^i^^'t') — 1 imposed to preserve 
the CMB signal is expressed in terms of spherical har- 
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monies as follows: 



-'oo 



= 0. (/ > 0) 



(10) 

(11) 



IV. DETERMINATION OF LINEAR WEIGHTS 

By minimizing the variance, we shall be able to derive 
equations that lead to the linear weights of minimum 



is real-valued, 
= (-1)"<™*- 



foreground. Since the function w (6 
^im obeys the reality condition wf_„ = 
Therefore, we need to determine only wfj-^-^ (m > 0). For 
computational convenience as well as to accommodate 
the reality condition, we define the real-valued spherical 
harmonic coefficients wf^^ as Re[wf^], Im[t(;f^] for m > 0, 
m < respectively. The constraints given by Eq. [10] and 



[TT]have the following forms for iwf^: 



j2wL = a>o). 



(12) 
(13) 



The linear weights of the minimum foreground will min- 
imize the variance X^la^ |o±2, lmP under the constraints 
Eg. I12l and ll3l The constrained minimization problem is 
solved conveniently via the Lagrange multiplier method 
[2l| . With the introduction of Lagrange multipliers 
Xim , it can be shown that the variance is minimized under 
the constraints Eq. [12] and [TS] when 



LM 



dwu' I 

I'm' 
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00- 



9w^'m' 



(14) 
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By using Eq.[8] it can be then shown that Eq.[T4]has the 
following form: 



E\k'k 



kin 



where af;'^,;„ is 



0, 



(15) 



(^^'m'lm — 2R-e 



Y^ll'il' .m\L,M)%{l,m,L,M) 



.LM 



(16) 



and 7fc(/i,mi,Z3,m3) is 

7fc(/l,mi,;3,TO3) + (-l)'"i7fe(?i,-TOi,;3,TO3) 

7fe(Zi,mi,;3,TO3) 
i [ikili, -miJs, TO3) - (-l)™i7fe(/i, mi, ^3, m^)] 



for mi > 0, mi = and mi < respectively, and 

7fc(/i,mi,/3,m3) = (17) 



E(-i)""' 



{2h + l){2h + l){2h + l) 



h h h 
mi m2 — m3 
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±2 T2 / "±2,'2m2 



Thus, the values of the linear weights of the minimum 
foreground can be found in terms of their Lagrange mul- 
tipliers Xi'm' by solving the system of simultaneous linear 
equations given by Eq. [15] The values of the Lagrange's 
multipliers Xi'm' can be easily determined by making the 
solutions of Eq. [15] satisfy the constraints Eq. [12] and [T3l 
Ea. ll2l[T5] and [T5l can be conveniently put in matrix form. 
For / frequency channels, in matrix notation Eq. 1151 be- 
comes 



A • w 



(18) 



where 



A ■/ 


— rv^'^ 

^I'm'lm 1 








= W^m, 




= Xi'm', 



for f = f{l'^ + l' + m') + fc', j = f{P + l + m) + k, and 
j" = l''^ + I' + m' . A is a n X n matrix, and w and L are 
column vectors of length n, and length n/ f respectively, 
where n is the total numbers of wf^,- 

In Eq. 1181 w is solved in terms of n/k undetermined 
Langrange multipliers, provided that A is invertible: 



w = A^n^L. 
The constraints given by Eq. [T^] and [T3] are: 
n • w = e. 



II is a J X n matrix, given by 



(19) 
(20) 



n, 



1 : /(*-l) + l<J </* 
: otherwise 



and e is a column vector of length n/ f, given by 



J = l 
J>1 







With Eq. and UM the undetermined n/f Langrange 
multipliers are given by: 



-(nA-in^)-ie, 



(21) 



Therefore, linear weights in spherical harmonic space is 
given as follows: 



A^n^(nA^n^)-ie. 



(22) 



It should be noted that Eq. [22] is not reduced to w = 
II~^e, since II and 11""" are not square matrices. 
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V. FOREGROUNDS REDUCTION FOR E AND 
B MODE POLARIZATION 

The polarization signal for the foregrounds as well as 
the CMB can be decomposed into E and B mode, where 
the E and B mode polarization maps at a frequency Vk 
are as follows: 

QE{vk,-n)±iUE{vk,-a) = - ^ a|_;,„ ±2rirn(n), 

Irn 

QB{vk,n)±iUB{vk,'^) = a| ±2ilm(n). 

Foreground sources for E and B mode polarizations can 
have distinct astrophysical origins. Therefore, frequency 
spectra for the foregrounds of E mode polarization might 
not be identical with those of the B mode polarization. 
Hence, we shall assume independent linear weights, 
and ui^, for the E and B mode, and find them respec- 
tively through the minimizations of <j%e ^'^'^ ^'bb^ where 

uIe = {\QE{h)±iUE{n)\^) 
ale = {\QB{n)±^UB{n)\') 
To determine w'^ ;„j and w% we replace i2m2 '^^ith 
-«l,i2m2 TiaB,i2m2' respectively, in Eq. [TT] After 
determing w'^ and w^, we then construct the following 
maps: 

e(n) = J2wUn) iQE{n,i^k) ± ^UEih,:yk)), (23) 

i 

/3(n) = ^w;|(n)(gB(n,i^fc)±zt/B(n,i.fc)). (24) 

i 

Using Eq. [?1 we may show 

^(n) = g^i"'^(n)±*C/r^fi) 

+ ^^«|(n) (Q|(n,i.fe)±zC/|(fi,j.fc)) 

i 
i 

/3(n) = g^i"'^(n)±iC/r^fi) 

i 
i 

Since the Hnear weights, w|;(n) and w^(n), are spatially 
varying functions, ^(n) contains the B mode polarization 
of the foregrounds and the noise, while /3(n) contains E 
mode polarization of foregrounds and noise. Therefore, 
we filter out the B mode polarization from ^(n) and the 
E mode polarization from /3(fi). We then reconstruct the 
CMB polarization map as follows: 

g(n)±z[/(n)=e(n)+/3(n), (25) 

where ^(n) and /3(n) are then filtered ^(n) and /3(n) re- 
spectively. 



VI. APPLICATION TO THE WMAP FIVE 
YEAR DATA 




0.0 30.0 micro K 



FIG. 1: the 1° FWHM smoothed HIL C5YR ma ps: tempera- 
ture + polarization, Q, U, and S = \/ Q'^ + (from top to 
bottom) 

We have applied our foreground reduction method to 
the WMAP five year polarization data Q. Since there 
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exists some anomalous excess power in the W band [3, 
[2^ . which is not fully understood, we do not use W band 
data in the reconstruction of the CMB polarization maps. 
We have considered only the E mode polarization and 
obtained E mode spherical harmonic coefficients of band 
maps as follows: 

j 

where a^"', are the E mode coefficients of the WMAP 
jth year single data in the fcth frequency channel, and 
the Bf are the beam transfer functions of the WMAP 
kth channel BUlil]. 

For the frequency channels of the multiple differencing 
assembly (Q and V band), we have used the average value 
of the multiple differencing assembly data. As discussed 
in Sec. IIIIl the cross terms in Eq. [6] bias the resulting 
CMB maps jl^]. As a part of HILC method, we had 
developed methods to reduce the cross term effect (see 
for details). However, in this work, we do not take 
special care of the cross term effect, since the cross term 
effect is negligible, compared to the noise contained in 
the resulting CMB map. 

When linear weights are obtained through variance 
minimization on noisy data, variance minimization pro- 
ceeds in the to minimize noise rather than foregrounds. 
Since the noise of the WMAP data dominates foreground 
signal on multipoles higher than ^ 60, we have used only 
afm {l ^ 60) in variance minimization (i.e. summation 
over I2 was done up to 60 in Eq. [T7]) . Though ideally the 
total number of wf^ may be as high as the total number 
of available af^ (i.e. the number of parameters may be 
as many as the number of data), we find that the number 
of w^^-^ , which keeps the matrices in Eq. [21] numerically 
non-singular, is much smaller than the ideal case (i.e. 
^cutoff <C 60). This may be attributed to the large band- 
width and the relatively small separation of the WMAP 
frequency channels, because in such configurations the 
frequency spectrum of foregrounds may not be numer- 
ically distinct enough over the channels. Starting with 
smallest number, we tried various Zcutoff and found the 
largest and stable Zcutoff for the WMAP polarization data 
is 5. Therefore, our linear weights Zcutoff = 5 is ineffective 
in cleaning foregrounds, whose frequency spectra vary on 
the angular scales smaller than ~ 36°. However, Fig. [1] 
shows that our CMB polarization maps do not contain 
significant foregrounds even inside the region correspond- 
ing to the WMAP team's Galactic cut. Besides that, we 
will be able to set ^cutoff as high as the coherence angu- 
lar scales of the Galactic foreground frequency spectra, 
when the polarization data from the upcoming Planck 
surveyor [2^ become available. 

As a consequence, linear weights of HILC maps are ob- 
tained through minimization of — X]l=2 m W±2,lm\'^ 
(see Eq.[S]). However, these linear weights can be applied 
to band maps containing higher multipoles. Hence, we 
used QE{n,i'k) ± illEin, v^) of multipoles (2 < / < 700) 



in the map reconstruction by Eq. [23] and [25] In Fig. [1] 
we show the CMB maps 'Harmonic Internal Linear Com- 
bination maps' (hereafter, HILC5YR). The top figure in 
Fig. [1] shows the CMB polarization as headless vectors 
whose length is proportional to polarization strength, 
while the underlying color-coded temperature map is 
from the HILC5YR temperature map [IjI- The second 
and the third figures show the Stokes parameter Q and 
U respectively, and the bottom figure shows the polariza- 
tion strength S = ^JQ^ + C/2. AU fi gures shown in Fig. [T] 
have been generated with 1° FWHM beam smoothing to 
suppress noise for illustrative purpose. 



200 




multipole I 

FIG. 2: TE correlation: solid smooth curve (the WMAP 
best-fit ACDM model), solid curve (HILC5YR), dotted curve 
(WMAP team's estimation) 



200 1 X WMAP 
150 + HILC 




-200 1 . . . 

50 100 150 200 250 
multipole I 



FIG. 3: X mark : the difference between the WMAP team's 
TE correlation and theory, -I- mark: the difference between 
HILC TE correlation and theory, solid curve (the error bars 
of the HILC TE estimation) , dashed curve (the error bars of 
the WMAP team's TE estimation) 
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In Fig. [2l we show the TE correlation of HILC maps 
and the estimation by the WMAP team. The theoretical 
TE correlation of the WMAP best-fit ACDM model is 
shown as a smooth solid curve in the same figure. Since 
noise in a temperature map and in a polarization map is 
uncorrelated, we have obtained temperature and E mode 
correlation by a simple quadratic estimator: 



21 



i-VRe 
+ 1 ^ 



(26) 



It can be shown that the variance of the estimator Cf^ 



\cF\') 



(27) 



21 + 



where C/ and Ni are, respectively, the true power spec- 
tra of the underlying CMB and noise distributions. In 
Fig. [21 we show the difference between the estimations 
and theory, and estimation errors. As shown in Fig. [3l 
HILC TE correlations are in agreement with theory and 
the WMAP team's estimations within expected errors. 
We have estimated Nf^ and N^^^ from the second term 
on the right hand side of Eq. IA4I and obtained the esti- 
mation errors by Eq. 1271 
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FIG. 4: E mode power spectrum: solid smooth curve(the 
WMAP best-fit ACDM model), solid curve(HILC5YR), dot- 
ted curve (WMAP team's power spectra) 

On the other hand, we cannot use a simple quadratic 
estimator for E mode power spectrum because of a noise 
bias. Hence we have derived a unbiased quadratic esti- 
mator of power spectra, which are similar to the WMAP 
team's cross power spectra, but in a convenient form 
to use with HILC method. The details on the unbi- 
ased quadratic estimator are given in Appendix [X] In 
Fig. m we show the E mode power spectrum obtained 
from HILC5YR with our estimator. We also show the 
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FIG. 5: X mark : the difference between the WMAP team's 
E mode power and theory, -I- mark: the difference between 
HILC E mode power and theory, solid curve (the error bars 
of the HILC E mode power estimation), dashed curve (the 
error bars of the WMAP team's E mode power estimation) 



theoretical E mode power spectrum of the WMAP best- 
fit ACDM model and the WMAP team's estimation ^ 
in Fig. [31 The large fiuctuation in comparison to the the- 
oretical prediction is attributed mostly to the estimation 
variance associated with noise (i.e. the third term on 
the right hand side of Eq. IA9p . Some estimated E mode 
power spectra are negative, because they are cross power 
spectra. In Fig. [3 we show the difference between the 
estimations and theory, and estimation errors. We have 
obtained the estimation errors of HILC E mode power 
by Eq. \M\ As shown in Fig. [3 we find that the HILC 
E mode power estimations are in agreement with the- 
ory and the WMAP team's estimations within expected 
errors. 

The linear weights of HILC5YR, which are continuous 
over the entire sky, are shown in Fig.[6l It is worth noting 
that the linear weights for Q and V bands are positive 
over the entire sky. We have computed the variance of our 



linear weights by wf = {21 + 1) ^ J^r 



k 12 
Im I 



to quantify 



the spatial variation of our linear weights on different 
angular scales. Fig. [3 shows that tends to decrease 
with increasing multipole with w§ being the highest. 



VII. COMPARISON WITH OTHER METHODS 

We restrict our discussion to the methods of currently 
available CMB polarization maps. Before proceeding to 
comparison, we would like to reiterate that the major 
goal of HILC approach is to reconstruct a whole-sky 
CMB map, which is important for the study of CMB 
polarization on large scales. 

The merits of HILC approach are that it naturally in- 
corporates spatial variability of the foreground frequency 
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FIG. 6: The HILC5YR linear weight for K, Ka, Q, and V 
band (from top to bottom) 




1 2 3 4 5 

multipole I 

FIG. 7: Variance of the HILC5YR linear weights 



spectrum in a natural way, and it barely relies on external 
information about foregrounds. Its major limitation is a 
poor performance for low Signal-to-Noise Ratio (SNR) 
data. 

The WMAP team produced foreground reduced maps 
by a template fitting method, whose model templates are 
derived from the K band of the WMAP data, and the 
Finkbeiner dust model 8 [2Q] . However, there are strong 
arguments against the interpretation of K band as syn- 
chrotron templates [2?, '2^ and it turns out that template 
models are not sufficient to make good fits simultaneously 
inside and outside the Galactic cut. Consequently, there 
is heavy foreground contamination inside the Galactic 
cut, which makes template- fitted maps unsuitable as a 
whole CMB sky map. 

The WMAP team also produced a low resolution CMB 
polarization map by the Markov Chain Monte Carlo 
(MCMC) method [9]. The merit of this approach is that 
it utilizes external information on polarized foreground. 
Its major limitations are that it does not work well within 
Galactic cut fij and its high computational cost makes it 
unable to produce a high resolution CMB map. 



VIII. DISCUSSION 

Using the WMAP 5 year polarization data, we have re- 
constructed whole-sky CMB polarization maps through 
Harmonic Linear Combination method. From our CMB 
polarization maps without any masking, we have ob- 
tained TE correlation and E mode power spectrum, 
which are consistent with the WMAP team's estimation. 

Because of the low Signal-to-Noise Ratio of the WMAP 
polarization data, our HILC polarization maps contain 
relatively significant noise. Besides that, the low Signal- 
to-Noise Ratio of the multi-frequency data seriously de- 
grades the foreground reduction of our method, which is 
relevant to all variants of ILC method. Therefore, we 
should warn that our CMB polarization maps should be 
used with some caution. However, the major aims of 
this work are to demonstrate HILC on polarization data 
in preparation for the future Planck data and to recon- 
struct whole-sky CMB polarizations maps, whose power 
spectra are similar to the WMAP team's estimation. 

The HILC method is not, in general, less effective in 
cleaning point sources than in the diffuse Galactic fore- 
grounds, because of assumed finite Zcutoff- However, we 
can effectively suppress the point-source contamination 
by applying a point-source filter derived by [1^ at the 
post-HILC stage. 

In comparison to other foreground cleaning method 
(e.g. template fitting methods), the effectiveness of our 
method improves sharply with an increase in Signal-to- 
Noise Ratio (SNR), the number of frequency channels 
and the angular resolution of the observation. Therefore, 
our method will allow us to reduce foregrounds effectively 
in the polarization data from the upcoming Planck sur- 
veyor [23, m. 
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The CMB polarization map, estimated power 
spectra and linear weights are available from 
http : //www . nbi . dk/ ~ jkim/hilc. 
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APPENDIX A: UNBIASED QUADRATIC 
ESTIMATOR OF POWER SPECTRA 



The estimator derived here is equally applicable to 
temperature, E mode and B mode power spectra. Since 
spherical harmonic coefficients of single- year maps con- 
sist of signal and noise terms: 



sL + N, 



(Al) 



the spherical harmonic coefficients of band maps, con- 
structed by averaging n single-year maps, are as follows: 
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where Sf^ denote the signal of fcth frequency channel 
and denotes the noise of the jth year data in the 

A:th frequency channel. Therefore, it is easy to show that 
the spherical harmonic coefficient a;„j of a HILC map is 
as follows: 



'S, 



Im 



(A2) 



where ® denotes spherical harmonic convolution given by 
Clebsch- Gordon relation. The spherical harmonic coeffi- 
cients of a HILC map may be then split into signal and 
noise terms: 



Sh 



3 

Im ' 



where Si^ = wf^ ® and N^^ = ® N^n^ ■ 

A simple quadratic estimator for power spectrum gives 

(2; + l)-i^|ai™p 

771 

\ m m j 

m j 

+i E E K.N''L + n'LnL 



m j>j' 



(A3) 



where the summation is over 1 < j' < n and j' < j < n. 
As seen from Eg. IA3[ the term l-ZV^-'^J^ is always positive, 
and hence makes Eq. IA3I biased. 

Now, consider the following quadratic estimator, which 
has an extra term: 



Ci = (2/ + i)-M^ 

\ m 



air, 



r^EEK 



1 2 



ii?{n — 1) 



(A4) 



m j>j' 



where 



aim = E ® ' = E ' 



k k 

Using Eq. EU [A3] and EH wc show that 
Ci = (2Z + l)-i \Y,\Sim\'' 

\ m 



m J 



'^^z^ E E NLN\^ + n^LnL ] 

m j>j' j 



(A5) 



Since the expectation value of Ci is 



lra\ / 7 



we find Ci unbiased. Now, let us consider the variance 
of the estimator Ci . We neglect residual foregrounds in 
signal Sim, assuming foreground reduction is very effec- 
tive. Since \Sim\'^/Ci follows distribution of 2Z-I-1 
degrees of freedom, the variance of the first term on the 
right hand side of Eq. IA5I is 



(lEl^'™l'l')=2(2/ + l)(Q)= 



(A6) 
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Taking into account the reality conditions Sim = 
(-l)'5r_„ and N}^ = (-l)'iV^7_^, we find the vari- 
ances of the second term and the third term on the right 
hand side of Eq. IA5I are respectively: 

(I J2 E SimN^L + SLKrf) = (2? + 1) • n . 4 . Q Ni, 

m j 

(A7) 

and 

m j>ji 

= {21 + 1)- ^ ^ ' -A - {Ni)\ (A8) 

where C; and Ni are the theoretical power spectrum of 
CMB and noise respectively. In deriving Eq. IA7I and 



IA8[ we have made the approximation {N^*i^Ni^,) oc 
Smm'j which is not true for real observations because of 
anisotropic noise. Hence our estimation on the variance 
of Ci will be slightly lower than the true values. Noting 
there is no correlation among the first, the second, and 
the third terms on the the right hand side of Eq. \A5\ we 
can show that the estimation variance of C'l is 

Zi + I n n(n — I) 

(A9) 

We can also obtain the same result with C; by applying 
the WMAP team's cross-power estimator to a pair of 
distinct single-year maps and averaging the estimation 
over all combinations. However, Ci by Eq. IA4I is more 
convenient for an HILC implementation. 
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